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Abstract 

This paper is devoted to the study of Hawking radiation as a tun- 
neling of charged fermions through event horizons of a pair of charged 
accelerating and rotating black holes with NUT parameter. We eval- 
uate tunneling probabilities of outgoing charged particles by using 
the semiclassical WKB approximation to the general covariant Dirac 
equation. The Hawking temperature corresponding to this pair of 
black holes is also investigated. For the zero NUT parameter, we find 
results consistent with those already available in the literature. 
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1 Introduction 



Hawking radiation [lj is a quantum mechanical process which can be de- 
scribed as follows: Firstly, a pair of positive and negative energy particles 
is created near the horizon. The negative energy particle falls inside the 
horizon, while the positive energy particle tunnels outside the horizon and 
rises in as Hawking radiation for the outside observer. Secondly, the pair 
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of particles is created just inside the horizon. The positive energy particles 
have ability to cross the energy barrier and tunnels to infinity, however the 
negative energy partner remains behind. This reduces the mass of black hole 
(BH) by losing energy in the form of Hawking radiation. The process in 
which particles have finite probability to cross the event horizon (that can- 
not be possible for classical particles) is called quantum tunneling. In this 
process, we evaluate tunneling probabilities that depend on the imaginary 
part of the action for emitted particles through the horizon. 

There are two methods to find the imaginary part of the action for the 
emitted particles. The first method was developed by Parikh and Wilczek [2] 
named as the null geodesic method by following the Kraus and Wilczek [3] 
technique. The second one is called Hamilton— Jacobi method which was first 
used by Angheben et al [I] for tunneling phenomenon to extend the complex 
path analysis of Srinivasan et al [5], [6] - The tunneling probability T for the 
emitted particle can be written as 

T oc exp[-2Im/], (1.1) 

where I is the classical action (h~ = 1). An approximation, called WKB 
approximation was used by Wentzel, Kramers and Brillouin (WKB) [7] to 
develop a relation between classical and quantum theory. This is used to 
solve the Schrodinger equation (h = leads to Hamilton— Jacobi equation 
in classical theory) in quantum mechanics. In both the above methods, the 
WKB approximation is used. 

The study of Hawking radiation as a tunneling phenomenon from various 
BHs has attracted many people. Page [8J calculated the emission rates from 
an uncharged, nonrotating BH for massless particles of spin ~, 1 and 2 by 
using the perturbation formalism. The same author j9] found numerical 
calculations of the emission rates of massless particles from a rotating BH. 
Kerner and Mann [10] extended the calculations of the tunneling process for 
the spin | particles emission by using the WKB approximation to the Dirac 
equation and found the tunneling probability from nonrotating BHs. Also, 
fermions tunneling is applied to a general nonrotating BH and recovered the 
corresponding Hawking temperature. The same authors [11] also discussed 
the tunneling of charged spin | fermions from a Kerr— Newman BH and 
obtained the Hawking temperature. 

In recent papers, the tunneling probabilities of incoming and outgoing 
scalar and charged/uncharged fermions from accelerating and rotating BHs 
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have been investigated [12]. Also, the thermodynamical properties of acceler- 
ating and rotating BHs with Newman— Unti— Tamburino (NUT) parameter 
have been studied [13] . Recently, we have also explored some work [T3] about 
the tunneling phenomenon for different BHs by using the above mentioned 
methods. In this paper, we extend this tunneling phenomenon of charged 
fermions from a pair of charged accelerating and rotating BHs with NUT 
parameter by using the procedure [10] . 

The paper is outlined as follows. In the next section, we briefly review 
Hamilton— Jacobi ansatz. Section 3 is devoted to explain the basic equations 
for a pair of accelerating and rotating BHs with NUT parameter having 
electric and magnetic charges. In section 4, we deal with the solution of 
the Dirac equation of charged particles in the background of these BHs and 
evaluate the tunneling probabilities as well as the corresponding temperature 
across the horizons. Section 5 provides an exact form of the action for the 
massless and massive particles. Finally, we summarize the results in the last 
section. 



2 Hamilton— Jacobi Ansatz: An Overview 

In this section, we briefly review Hawking radiation as tunneling through 
Hamilton— Jacobi ansatz [3l[3],[T0]. The metric for a general (non-extremal) 
BH is 

ds 1 = -f(r)dt 2 H — — — - + Cirthadx'dx 3 . (2.1) 

g{r) 

Assume a scalar particle moving in this BH. Since we are interested only 
in the leading terms within the semiclassical approximation, we neglect the 
effects of the particle self-gravitation. Thus the quantity which satisfies the 
relativistic Hamilton— Jacobi equation is the classical action I 

g^d^WJ + m 2 = (2.2) 

which implies that 

" + 9(r)(d r I) 2 + -^rfiW + m 2 = 0. (2.3) 

Due to the symmetries of the metric, there exists a solution of the form 

1 = -Et + W{r) + J(x { ), (2.4) 
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where d t I = —E, d r I = W'(r), dil = Jj and the Jj's are constants. E 
corresponds to energy as d t is the timelike Killing vector. Solving for W(r), 
it follows that 



W ' ±(rH± /77feMf 2 - /(r) ( ra2 + Z ff)' (2 ' 5) 

Here, W + and W_ correspond to scalar particles moving away (outgoing) and 
moving towards the BH (incoming), orespectively. Notice that the imaginary 
part of the action can only be due to the pole at the horizon. 
The probabilities of crossing the horizon each way are 

2 2 
Probfout] oc exp[— — Imil = exp[— — ImW + ], (2.6) 
a a 

2 2 

Probfin] oc exp[— — Im/] = exp[— — ImW_]. (2-7) 
n n 

Since W + = —W-, we obtain the probability of a particle tunneling from 
inside to outside the horizon 

4 

r cxexp[--ImW r +]. (2.8) 

We take H to be unity and also drop the '+' subscript from W. Integrating 
around the pole at the horizon leads to the result [15] 

W= - . (2.9) 

Consequently, the tunneling probability becomes 

47T 

T = exp[ E] (2.10) 

which leads to the usual Hawking temperature 

T H = . (2.11) 
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3 Accelerating and Rotating NUT Solutions 

Black holes are the most important predictions of general relativity. The 
discovery of Schwarzschild BH is followed by the extension to the electrically 
charged version, Reissner— Nordstrom (RN) BH. This is followed by the dis- 
covery of rotating versions, Kerr and Kerr— Newman BHs [16] . The research 
in the BH area has further been extended by adding different sources like 
cosmological constant, an acceleration as well as a NUT parameter. Black 
hole solutions with these extensions belong to type-D class [PT] . 

In general, the NUT parameter is associated with the gravitomagnetic 
monopole parameter of the central mass, or a twisting property of the sur- 
rounding spacetime but its exact physical meaning could not be ascertained. 
Recently, higher dimensional generalization of the Kerr— NUT— (anti) de Sit- 
ter spacetime [18] and its physical significance [19] is investigated. For this 
BH, the dominance of the NUT parameter over the rotation parameter leaves 
the spacetime free of curvature singularities and the corresponding solution is 
named as NUT-like solution. However, if the rotation parameter dominates 
the NUT parameter, the solution is Kerr-like and a ring curvature singularity 
forms. This kind of behavior on the singularity structure is independent of 
the presence of the cosmological constant. 

There are many BH solutions which incorporate the NUT parameter and 
investigate its physical effect in the space of colliding waves. Exact interpre- 
tation of the NUT parameter becomes possible when a static Schwarzschild 
mass is immersed in a stationary, source free electromagnetic universe [20] . In 
this case, the NUT parameter is related to the twist of the electromagnetic 
universe excluding the central Schwarzschild mass. In the absence of this 
field, it reduces to the twist of the vacuum spacetime. Thus, the twist of the 
environmental space couples with the mass of the source to generate NUT 
parameter. It would be interesting to explore the tunneling phenomenon 
from such a BH. 

The family of solutions which represents a pair of charged accelerating 
and rotating BHs with a nonzero NUT parameter can be written as [21] 




(3.1) 



5 



where 



n 
Q 

P 



a 



1 (1 + a cos (9)r, p = r 2 + (I + a cos 0) , 

2;. .2 ..2,/, , .-,«".., o 1 r.. ~2 



[cu 2 k + e 2 + g 2 )(l + 2—r) - 2Mr + 



a(a — I) 
1 + — i -r 



■jj 



UJ 

a (a + I) 

j 

UJ 



I 2 



sin 2 ^(1 — a 3 cos 9 — cos 2 - 
a 2 al 



Psin 2 ^, 



2 — M-4 

a; 



(w^fc + e z + 



a 4 



a 2 a 2 



UJ Z 



uj 2 k + e 2 + g 2 ). 



Here M represents a pair of BHs mass, e and g indicate electric and magnetic 
charges, respectively, while a shows a BH rotation and I is a NUT parameter 
[22] ■ The continuous parameter a represents acceleration of the sources. 
Also, the rotation parameter uj is related to the rotation of the sources and 
k is given by 



uj- 



3n 2 / 2 )/-' 1 -2-M-3^(e 2 + ^ 2 ), 



UJ 



UJ^ 



(3.2) 



where a, uj, M, e, g and k are arbitrary real parameters. 

Notice that this type of BH involves acceleration a while twisting be- 
havior of the sources is proportional to the parameter uj which is related to 
both the Kerr-like rotation parameter a and the NUT parameter I. Thus, 
a, M, e, g, I, a vary independently while uj depends on nonzero value of rota- 
tion parameters I or a. For a = 0, this reduces to the Kerr— Newman— NUT 
solution and Eq. fl3.2p will become uj 2 k = a 2 — I 2 . In the absence of NUT 
parameter I = 0, it reduces to the pair of charged accelerating and rotating 
BHs. Further, a = leads to the Kerr— Newman BH and a = yields the RN 
BH. In addition, if e = = g, we have a Schwarzschild BH while I = = a 
leads to the C-metric. 

The metric (13.1 p can also be written in a more suitable form as 



ds 2 = -/(r,^)dt 2 + 



dr 2 
9(r,9) 



+ S(r, #)d# 2 + K{r, 9)d(j) 2 - 2H(r, 0)dtd<j>, (3.3) 
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where f(r,9), g(r,9), S(r, 9), K(r,9) and H(r,9) can be denned as follows: 



f(r,9) 
K(r, 9) 
H(r,6) 



Q - Pa 2 sin 2 9 

phT 2 



g(r,9) = ^, E(r,0) 
P 2 



p 2 Q? 
1 

^fi 2 



n 2 p 

9 



sin 2 9P[r 2 + (a + /) 2 ] 2 - Q(a sin 2 9 + 41 sin 2 -) 2 



sin 2 9Pa[r 2 + (a + I) 2 ] - Q(a sin 2 9 + 41 sin 1 



9, 



(3.4) 
(3.5) 
.(3.6) 



The electromagnetic vector potential for these BHs is given by [2 

1 



.4 



-[-er[adt - d</>{(Z + a) 2 - (I 2 + a 2 cos 2 9 



a[r 2 + (I + acos9) 2 } 
+ 2lacos9)}] - g(l + a cos 9)[adt-d(j){r 2 + (I + a) 2 }]]. (3.7) 

The event horizons are obtained for g(r, 9) = = [IT] , where A(r) = 



^. This implies that A(r) = = Q, yielding the horizon radii 



CO 



E(r,9) 

e b 
a 2 -/ 2 



ai 



a(a + /) 



a(a — /) 



uj 2 k 



[(u 2 k + e 2 + g 2 )- 

UJ 



-M] ± \l [(u 2 k + e 2 + g 2 )—- M} 2 - -^L^(u 2 k + e 2 + g 2 ^ 



UJ 



I 2 



(3.8) 



where r± represent the outer and inner horizons, respectively, such that 
[(uj 2 k + e 2 + g 2 )f - M] 2 - ^f%(oj 2 k + e 2 + g 2 ) > 0, r ai and r a2 are ac- 
celeration horizons. The angular velocity at the BH horizon can be defined 

as 

H(r+,9) a 



n 



H 



The inverse function of f(r, 9) is 



F(r,9) = f(r,9) + 



r\ + {a + l) 2 ' 
H 2 (r,9) 



K(r, 9) ' 

For these BHs, this can be written in the following form: 

PQ sin 2 9p 2 



(3.9) 



(3.10) 



F(r,9) 



n 2 [sin 2 9P[r 2 + (a + I) 2 ] 2 - Q(a sin 2 9 + 41 sin 2 §) 2 ] 



(3.11) 



In terms of A(r) and E, we can write the inverse function of f(r, 9) as 

P 2 A(r)£(r,6>) 



F(r, 



[r 2 + (a + /) 2 ] 2 - A(r) sin 2 9[a + T 



(3.12) 



+COS ( 
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4 Tunneling of Charged Fermions 

In order to study charged fermions tunneling of mass m from a pair of accel- 
erating and rotating BHs with NUT parameter having electric and magnetic 
charges, the covariant Dirac equation with electric charge q is given by [2] 

t-f (D„ - jA^ * + jr* = 0, fi = 0, 1, 2, 3 

where A„ is the 4-potential, \& is the wave function and 



(4.1) 



D, 



(4.2) 



The antisymmetric property of the Dirac matrices [12], i.e., [7", 7^] = for 
a = (3 and [7°, 7^] = — [7 /3 ,7 a ] for a ^ (3, reduces the Dirac equation (14. ip 
in the following form: 



*y a 



"7 



ft 

The Dirac matrices are given by 



m 



I H' • — M' 0. 



(4.3) 



7 



7 



7 V 



1 
1 



1 
-1 

a 1 
a 1 

o 2 
a 2 



7' 



a 3 
a 3 



+ 



H(r,> 



/ F(r,8)K(r,6 

where a 1 (i = 1, 2, 3) are the Pauli sigma matrices defined as 



t 
-t 



a 



1 

1 



a 



-t 
t 



a 



1 
-1 



(4.4) 



The spinor wave function \I/ (related to the particle's action) has two spin 
states: spin-up (in +ve r-direction) and spin-down (in -ve r-direction) . For 
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the spin-up and spin-down particle's solution, we assume [TO] 



A(t,r,9,< 


B{t,r,0,, 




C{t,r,9,> 


D(t,r,9,. 



cxp 



exp 



(4.5) 



(4.6) 



where I-jyj, denote the emitted spin- up/spin-down particle's action, respec- 
tively. Note that we shall only analyze the spin-up case as the spin-down 
case is just analogous. 

The particle's action is given by the following ansatz (12. 4ft : 



I t = -Et + J<t> + W(r,0), 



(4.7) 



where E, J and W are energy, angular momentum and arbitrary function, 
respectively. Using this ansatz into the Dirac equation with iA = B, lB = A 
and Taylor's expansion of F(r, 9) near the event horizon, it follows that 



B 



B 



-E + Q H J 



qer 



■s/(r - r + )d r F(r+ 



+ V(r-r + )d r g(r + ,e)(d r W) 



+ niA = 0, 
(4.8) 



'Pfi 2 (r H 



P 2 (r + ,0) 



+ 



x 



id e W) 



sin 2 9P[r 2 + (a + I) 2 } 2 - Q(a sin 2 9 + 41 sin 2 |) 2 

\er [{I + a) 2 - (I 2 + a 2 cos 2 9 + 2la cos 9)} 
~ q [ a[r 2 + (/ + a cos 9) 2 } 

g(l + acos#)[r 2 + (I + a) 2 ] 



a[r 2 + (I + a cos ( 



0. 



(4.9) 
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.4 



A 



-E + Q H J + 



qer 



y/{r-r + )d r F{r + ,9) 



- V(r-r + )d r g(r + ,0)(d r W) 



mB = 0, 
(4.10) 



(d W) 



+ 



L P (r + ,B)Q(r + ,e) 



sin 2 9P[r 2 + (a + I) 2 ] 2 - Q(a sin 2 9 + 41 sin 2 § ) 2 
er[{l + a) 2 - {I 2 + a 2 cos 2 9 + 2la cos 9)] 



a[r 2 + (1 + a cos #) 2 ] 
#(7 + acos6>)[r 2 + (/ + a) 2 ] ^ _ 



x { J-q 

^ a[r 2 + (I + a cos 9) 2 
The arbitrary function W(r, 9) can be separated as follows: [TT 

W(r,0) = R(r) + Q(9). 



(4.11) 



(4.12) 



First we solve Eqs. ( 14.8(1 — ( 14.111) for the massless case (m — 0). Using the 
above separation, Eqs. (14. 8 j) and (14. lOf) reduce to 



-E + Q H J + 



qer 



- B 



A 



[ri+(*4-Q a ] 



y/{r- r + )d r F(r + , 



qer 



+ y / (r-r + )d r g(r + ,9)R'(r) 



E + Q H J + [r 2 +(a+0 2] n — — — 

V(r-r +)dr F(r + ,9) ~ ^-^'^ W 



0, (4.13) 
= 0, (4.14) 



implying that 



R'(r) = R'Ar) 



-R' (r) 



[r\ + (a + I) 2 ] 



(r - r + )[l + 



a(a—l) 



r+][l- 



q(a+Q 



ger + 



X 



2[^( a ;2 A; + e 2 + ^2 ) _ M+ _Q_ r+] 



(4.15) 



where i? + and i?_ correspond to the outgoing and incoming solutions, re- 
spectively. This equation represents the pole at the horizon, r = r + . 

In order to calculate the Hawking temperature by using the tunneling 
approach, we need to regularize the singularity by specifying a suitable com- 
plex contour to bypass the pole. In standard coordinate representation, for 
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outgoing particles (form inside of the horizon to outside), we should take the 
contour to be an infinitesimal semicircle below the pole r = r + . Similarly, for 
the ingoing particles (from outside to inside), the contour is above the pole. 
Here, for the purpose of calculating the semiclassical tunneling probability, 
we need to multiply the resulting wave equation by its complex conjugate. 
In this way, the part of trajectory that starts from outside of the BH and 
continues to the observer, will not contribute to the calculation of the fi- 
nal tunneling probability and can be ignored (since it will be entirely real). 
Therefore, the only part of the wave equation (trajectory) that contributes 
to the tunneling probability is the contour around the BH horizon. Instead 
we choose a mathematically equivalent convention that the outgoing contour 
is in the lower half plane and so do not multiply by a negative sign [11] . 
Integrating Eq.f l4.15p around the pole, we obtain 



Rjr) 



-i?_(r) 



7ri[rl + (a + If 



[l + ^ + ][l-^r + ] 

[E-n H j- [r 2 q + e [ a + +l) 2 ]S 



x 



2 [sL(u*k + t* + g*) 
The imaginary parts of R+ and R- yield 



M + 



. (4.16) 



ImiL 



-ImR. 



77[rl + (a + /) 2 ] 



[1 + 



a(a— I) 



r+][l 



a(a+l) 



[e-q h j- } 



r 2 + + (a+l) 2 ]i 



2 [^ k + e 2 + g 2 ) _ M+ _gk u 



(4.17) 



Thus the particle's tunneling probability from inside to outside the horizon 



is 



Prob [out] 
Prob [in] 



exp[-2(Imi?+ + Im0)] 



exp[-2(Im J R_ + Im0)] 

-2n[rl + (a + l) 2 ] 



exp[— 4Imi? H 



exp 



r-i 



_|_ a(a-l) r j ^ _ a(a+l) 

[E-n H j- qer+UJ 1 



X 



[r 2 +(a+«) 2 ]J 



[ ^ u 2 k + e 2 +g 2 ) _ M+ ^_ r+] 



(4.18) 
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Here, T is the tunneling probability for the classically forbidden trajec- 
tories of the s-waves coming from inside to outside the horizon. Using the 
WKB approximation, T is given by Eq.(l.l) in terms of classical action / of 
Dirac particles tunneling across the BH horizon as trajectories up to leading 
order in h. Thus, for calculating the Hawking temperature, we expand the 
action in terms of particles energy E, i.e., 21 = j3E + 0(E 2 ) so that the 
Hawking temperature is recovered at linear order given by 

T ~ exp[-2i] ~ exp (4.19) 

This shows that the emission rate in the tunneling approach, up to first order 
in E, retrieves the Boltzmann factor, exp[— (3E], where /3 — ^- [15]. The 
higher-order terms represent the self-interaction effects resulting from the 
energy conservation. 

Semiclassically, for the charged rotating BH solution, the energy E of the 
emitted particle in the Boltzmann factor should be replaced by E — VLhJ due 
to the presence of the ergosphere. The corresponding Hawking temperature 
at the event horizon takes the form 

r [^k + e 2 + fl 2 ) - M + ^r + ] [1 + ggrj [1 - ^±V] ' 
H ~ [ 2vr[r 2 h + (a + /) 2 ] 

(4.20) 

Thus, the Hawking temperature of Dirac particles tunneling from the event 
horizon of the charged accelerating and rotating BHs with NUT parameter 
is well described via fermions tunneling method. When 1 = 0, k = 1 in 
Eq. ( 14.20)) . the Hawking temperature of the accelerating and rotating BHs 
with electric and magnetic charges is recovered [12] . When a = 0, the 
Hawking temperature of the charged accelerating and rotating BHs is reduced 
to the temperature of non-accelerating BHs [13J. For 1 = 0, k = 1, a = 
in Eq. (14.201) . the Hawking temperature of the Kerr— Newman BH [TTJ is 
obtained which is further reduced to the temperature of the RN BH (for 
a = 0). Finally, in the absence of charge, the temperature exactly reduces to 
the Hawking temperature of the Schwarzschild BH |25j . 

For the massive case (m ^ 0), following the same steps, we can obtain 
the same temperature (14.201) . Thus the behavior of massive particles near 
the BH horizon is the same as that for the massless particles. For / = 0, the 
tunneling probability (14.181) reduces to the form for the pair of accelerating 
and rotating BHs [12]. For a = 0, we recover the tunneling probability for 
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the Kerr— Newman BH [IT]. In the absence of rotation, this yields the same 
result as for the RN BH [21]. 

Now we explore the tunneling probability of charged massive and mass- 
less fermions from the acceleration horizon r ai given in Eq.f l3.8p . The cor- 
responding set of Eqs. fl4.8p — (14. lip for the outgoing and incoming fermions, 
respectively, yield 



R + (r) 



-RJr) 



7n[r 2 ai + (a + l) 2 } 
[I + - P)r ai ] 
[E — Q a J ^ 



x 



^!+(a+0 2 ] J 



2^ + e* + 02)(l + - 2Mr ai + f^r\ x \ 



.(4.21) 



The tunneling probability can be written as 

-2n[r 2 ai + (a + lf\ 



exp 



Z[l + Z(a* - P)r ai ] 

[E - Vl a J 



qer ai 



X 



[rl^+l) 2 }"- 



[ {uJ 2 k + e 2 + g 2 ){1 + 2S^=L) _ 2Mr ai + - a 



2_p 



(4.22) 



Consequently, the Hawking temperature for the acceleration horizon turns 
out to be 



T, 



^ [{u 2 k + e 2 + g 2 ){1 + 



2alr a 



2M S + ^|][I + ^ 2 -I 2 )r Ql ] 



Mr 2 ai + (a + If 



(4.23) 



5 Action for the Emitted Particles 

In order to obtain the explicit expression for the action 1^ in the spin-up 
case, we solve Eqs. fl4.8p — (14. lip near the BH horizon. Using Eq. fl4.12p . we 
can write Eq.f l4.8p in the form 

i?(r) , mA - E + y + rf&m (51) 

By/(r -r + )d r g(r + ,9) (r - r + )^d r F(r + ,9)d r g(r+,9) 
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Integrating with respect to r, we get 



mA 



R(r) = R + (r) = [ 



-.dr 



[H+(a+0 2 ] 



y/d r F(r + ,6)d r g(r + ,9) 
Similarly, for the incoming particles, Eq.f l4.10l) yields 

mB 



ln(r — r + ). 



R(r) = R_(r) = J 



A^J(r - r+)d r g(r 
(~E + tt H J + 



zdr 



qerj, 



y/d r F(r + ,6)d r g(r + ,0) 
Using Eq.glSD, Eqs.gSJ) and (ICTj) imply that 



ln(r — r + ). 



+ 



P 2 ( r +> 9 ) ^/sin 2 6P[r\ + (a + I) 2 ] 2 

j er + [(I + a) 2 - (I 2 + a 2 cos 2 9 + 11a cos i 
J ~ q \ a[r\ + (/ + a cos 9) 2 } 

g(l + a cos 9)[r\ + (I + a) 2 ] v 



(5.2) 



(5.3) 



(5.4) 



a[r 2 + + (I + acos9) 2 ] 
Substituting the values of p and P, after some manipulation, it follows 
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that 



+ 



lcl 2 sin 9 J + iqer + a sin 9 



[l -2cos9 ^a^-M 



[r 2 +(a + /) 2 ] 

- 2a 2 l — (u 2 k + e 2 + g 2 ))- cos 2 \-a 2 — (co 2 k + e 2 + g 2 ) 

CO 2 J { co 2 

— taJ + iqqil + a cos 60 r „ . r a , , 

+ -\l-2cos9\a-M 

a sin 9 L La; 

- 2a 2 /^-(^ + e 2 + g 2 )} - cos 2 {-a 2 — (u 2 k + e 2 + g 2 ) 

co 2 ) i co 2 

2lat(l — cos 9) J + iqer + 2l(l — cos 9) 



I - 2cos9 \a-M 

CO 



{■ 



sin 9[r% + (a + I) 2 ] 
2e?l—Aco 2 k + e 2 + g 2 ))- cos 2 {-a 2 —Aco 2 k + e 2 + g 2 ) 

L0 2 J \ CO 2 



Integrating with respect to 9, we have 

4gaer + + Ja 2 ] 2%J + ger + ] 

[r 2 +(a + /) 2 ] 1+ 2 + [r 2 +(a + /) 2 ] 3 ' 

where Ji, J 2 and I 3 are given as follows 

h = [sin9 \l -2cos9la-M -2a 2 lA,(oo 2 k + e 2 + g 2 )} 
J i I to co 2 ) 



, a 



- cos 2 9 {-a 2 ^(u 2 k + e 2 + g 2 ) 



CO' 



-1 



d9, 



h. = 



L {<19(1 + a cos 9) — a J} 
asin6> 



1 - 2cos# 



la—M 
I co 



- 2a 2 l—(u 2 k + e 2 + g 2 )) - cos 2 {-a 2 — (uj 2 k + e 2 + g 2 )) 

CO 2 ) i CO 2 J 



/ 



1 — cos 9 
sin 9 



1 - 2 cos 9 I a-M - 2a 2 l^(co 2 k + e 2 + g 2 )\ 

I CO CO 2 ) 



- cos 2 9 < -a 2 — (co 2 k + e 2 + g 2 



CO 



-I -1 



c\9. 
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Solving these integrals, we obtain after some algebra 



h = 



x In 



2a^{M - 2a^(u 2 k + e 2 + g 2 )f - (\/u 2 k + e 2 + g*) 2 _ 
1 - a- cosO I M - 2a-(u 2 k + e 2 + g 2 ) 



+ \J^M- 2a l -(u 2 k + e 2 + # 2 ) J - (v^A 7 + e 2 + # 2 ) 2 j j 
x (l - a— cos 6 <M - 2a— (u 2 k + e 2 + g 2 ) 
- - 2a l -(u 2 k + e 2 + g 2 ) J - ( v^A; + e 2 + # 2 ) 2 



(5.10) 



L x In 



1 — a— cos # < M — 2a-(u 2 k + e 2 + g A 

CO I U) 



+ y / |M-2a-J-(w 2 A; + e 2 + ^ 2 )| - ( v / ^ : + e 2 + ^ 2 ) 2 ► j 
x ( 1 - a-cos# { M -2a-(u 2 k + e 2 + g 2 ) 

- ]JJm - 2a l -(u 2 k + e 2 + # 2 )j - ( v^A; + e 2 + # 2 ) 2 
1 - 2 cos ( a-M - 2a 2 /4r(a; 2 A; + e 2 + c/ 2 ) ) 



+ L 2 ln 



, a 



+ a 2 — cos 2 6(u 2 k + e 2 + g 2 ) 



+ L 3 ln[l - cos 9] 



+ L 4 ln[l + cos0], 



(5.11) 
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where 

u = 



x 



+ 



X 



| 2a ^ { M - ^-{yj^k + e 2 + # 2 ) j - (y/u 2 k + e 2 +^| 



1 + a 2 — (w 2 £; + e 2 + # 2 ) -4a 



2 



x M-2a-(w 2 £; + e 2 + # 2 



M 



2a 2 l—(u 2 k + e 2 + 5 2 )V - a^—(u 2 k + e 2 + <? 2 ) 2 
w 2 / a> 4 

5 2 ) ? + ^ i-a-M + 2a 2 l^(u 2 k + e 2 + a 2 ) 
J I a; ur 



a —zioj k + e 

47 4 

at a 



2a 2 / a 2 



2\2 



+ -(u; 2 £; + e 2 + (? 2 ) 



a or 
a 2 / a 2 



M-2a-(u; 2 A; + e 2 + £ 2 ) 



Au 2 k + e 2 + g 2 ) + a A — (u 2 k + e 2 + g' 
a u A ur 



M -2a-(u 2 k + e 2 + g 2 



(5.12) 
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Lo = 



l + a 2 -Ju 2 k + e 2 + g 2 )} 2 -4a 2 - 2 



OJ' 



x ^M-2a l -{u 2 k + e 2 + g 2 )^ 

x \ L j\ a -M-2a 2 l^(u 2 k + e 2 + g 2 )} 
^( 2 i (a*M-2a 2 l^Ju 2 k + e 2 + g 2 )) 



, a 



+ l + a 2 ^(w 2 A; + e 2 + # 2 ) 



^ ( - + 1 ) - i>J 



[l -2|a^M-2a 2 / 



x ^(w 2 /c + e 2 + /))+a 2 ^(w 2 /c + e 2 + # 2 ) 



L 4 = - t J + tggr 1 

2 [ V a 

+ 2 {a-M-2a 2 l—(u 2 k + e 2 + g 2 )) 
La; a; 2 J 

and J 3 can be obtained as 

I 



l + a 2 -(u 2 k + e 2 + g 2 ) 



-i 



iViln 



1 — a— cosO 1m — 2a-(u 2 k + e 2 + g 2 ) 
uj I a; 



+ 



|m - 2a l -(u 2 k + e 2 + # 2 )J - ( + e 2 + # 2 ) 2 



x ( 1 — a— cos # < M — 2a—(u 2 k + e 2 + g' 



+ N 2 \n 



M - 2a-{u 2 k + e 2 + # 2 )J - (^u 2 k + e 2 + g 2 ) 2 

1 - 2 cos 6 {a-M - 2a 2 l—(u 2 k + e 2 + g 2 ) } 
I u a; 2 J 



+ a 2 -^cos 2 6»(a; 2 A; + e 2 + g 2 ) 



+ N 3 \n[l + cos 9], 



where 



a±M - 2la 2 ^(u 2 k + e 2 + g 2 ) + a 2 ^{u 2 k + e 2 + g 2 ) 



1 



2a^sJ {M - 2a l -{yj 2 k + e 2 + g 2 )} 2 - {^u 2 k + e 2 + g^) 2 

2 

1 + a 2 ^(u 2 k + e 2 + g 2 ) + 2a{-M 



X 



N 2 




(5.17) 



iV 3 



2a 2 l-Ju 2 + e 2 + g 2 )) 
a 2 

- 1 + a 2 — (u 2 k + e 2 - 



CO 




(5.18) 



} 



i 



2a 2 l—Ju 2 k + e 2 + g 2 ) 



(5.19) 



Inserting these values of the integrals in Eq. fl5.6p and using Eq.f l4.12p . we 
obtain W(r, 6). This leads one to evaluate the action for the outgoing massive 
particles. For m = 0, this reduces to the action for the massless particles. 
Similarly, we can determine the action for the incoming massive and massless 
particles. 



In this paper, we have used semiclassical WKB approximation to study tun- 
neling of charged fermions from a pair of accelerating and rotating BHs hav- 
ing electric and magnetic charges, together with a NUT parameter. Clas- 
sically, a particle can only fall inside the horizon while in the semiclassical 
approach, the horizon plays a role of two way energy barrier for a pair of 
positive and negative energy particles. The positive energy particles have 
ability to tunnel outside the event horizon, which contradicts classical ap- 
proach. Thus, we have considered tunneling probabilities for both incoming 
as well as outgoing particles. Relating these tunneling probabilities with the 
Boltzmann factor exp[— f3E] for emission at the Hawking temperature, we 
can recover the corresponding Hawking temperature for this pair of BHs at 
event horizons. 



6 Outlook 
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The tunneling probabilities of outgoing/incoming charged fermions do not 
depend upon mass of the fermions but only its charge. The corresponding 
Hawking temperature depends upon mass, acceleration and rotation param- 
eters and also NUT parameter as well as electric and magnetic charges of the 
pair of BHs. Equations for the spin-down case are of the same form as for the 
spin-up case except for a negative sign. For both massive and massless cases, 
the Hawking temperature implies that both spin-up and spin-down particles 
are emitted at the same rate [TO] . 

It is worth mentioning here that in the absence of the NUT parameter I = 
0, all the results reduce to the results of the pair of charged accelerating and 
rotating BHs [T2]. Further, a = provides the results of the Kerr— Newman 
BH [11J and a = yields the results of the RN BH [21]. 
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